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Executive  Summary 


Background. 

The  application  of  cryogenic  refrigerators  and  similar  cryogenic  heat  sinks  to  the  cooling 
of  scanning  detectors  on  three  axis  stabilized  spacecraft  can  be  most  effectively 
implemented  by  mounting  the  cryogenic  heat  sink  on  the  despun  portion  of  the  vehicle 
and  cooling  the  scanning  sensor  via  a  rotary  cryogenic  thermal  coupling. 

Phase  I  Effort. 

The  principles  and  practices  required  for  the  design  and  fabrication  of  a  Rotary  Cryogenic 
Thermal  Coupling  were  developed  and  demonstrated. 

A  detailed  mathematical  analysis  of  the  radiant  heat  transfer  occurring  within  the  basic 
Rotary  Cryogenic  Thermal  Coupling  module  was  prepared,  including  a  mapping  of  the 
module  disk  fin  effectiveness  versus  the  disk  radius  ratio  and  a  component  dimensionless 
characterization  parameter,  X,  for  a  specific  temperature  regime. 

A  computer  program  was  written  to  parametrically  size  a  Rotary  Thermal  Coupler  based 
on  the  mathematical  analysis.  A  printout  of  the  results  of  several  different  cases  was 
presented. 

The  thermal  analysis  demonstrated  that  the  thermal  heat  leak  parasitic  losses  can  be  kept 
within  a  reasonable  percentage  of  the  transferred  heat  by  incorporating  active  thermal 
shielding.  The  technology  is  applicable  down  to  10  K. 

The  preliminary  design  analysis  concluded  that  the  unproven,  critical  component  portion 
of  the  Rotary  Cryogenic  Thermal  Coupling  was  the  basic  radiation  coupled  disk 
assembly.  Several  alternate  materials  and  processes  approaches  were  considered.  It  was 
concluded  that  a  concept  based  on  vacuum  brazed  aluminum  fabrication  was  the  best 
option.  A  proof  of  concept  module  was  designed  and  fabricated.  This  task  demonstrated 
that  the  module  components  can  be  fabricated. 

Phase  II. 


A  proposal  for  Phase  II  was  not  submitted  because  there  was  no  formal  request  issued  by 
the  contracting  agency. 
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Rotary  Cryogenic  Thermal  Coupling 

1 .0  Identification  and  Significance  of  the  Problem 

Conceptually,  the  application  of  cryogenic  refrigerators  and  similar  cryogenic  heat  sinks 
to  the  cooling  of  scanning  detectors  on  three  axis  stabilized  spacecraft  can  be  most 
effectively  implemented  by  mounting  the  cryogenic  heat  sink  on  the  despun  portion  of 
the  vehicle  and  cooling  the  scanning  sensor  via  a  rotary  cryogenic  thermal  coupling.  If 
the  cryogenic  heat  sink  is  to  be  either  a  cryogenic  radiator  or  an  expendable  cryogenic 
fluid,  vehicle  geometry  constraints  usually  show  that  it  is  impractical  to  mount  the  heat 
sink  device  on  the  scanning  sensor  segment.  If  a  cryocooler  is  to  be  used,  installing  the 
cryogenic  refrigerator  (and  its  heat  rejection  thermal  radiator)  on  the  despun  segment  of 
the  vehicle,  which  could  be  accomplished  if  a  rotary  cryogenic  thermal  coupler  were 
utilized,  would  result  in  a  reduction  in  weight,  power,  and  sensor  vibration  and  an 
increase  in  reliability. 

Unfortunately,  rotary  cryogenic  thermal  coupling  technology  has  yet  to  be 
demonstrated.  This  Phase  I  SBIR  research  and  development  activity  establishes  a 
basis  for  the  design  of  such  devices. 

1.1  Background 

Consideration  of  the  general  problem  of  totally  integrated  spacecraft  borne  cryogenic 
thermal  management  systems  indicates  that  rotary  thermal  couplings  are  required  for 
some  system  implementations.  A  specific  example  would  be  to  provide  for  the  transfer 
of  heat  from  a  rotating  (scanning)  focal  plane  to  a  cryogenic  heat  sink  located  on  de¬ 
spun  segment  of  the  spacecraft.  The  basic  requirements  for  an  acceptable  rotary 
thermal  coupling  would  be: 

1.  low  thermal  impedance; 

2.  negligible  internal  heat  generation; 
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3.  relatively  low  parasitic  heat  leak; 

4.  long  service  life; 

5.  high  reliability;  ■ 

6.  relatively  low  mass; 

7.  low  mechanical  noise; 

8.  robust  design  to  withstand  ground  handling  shock  and  launch  vibration; 

9.  convenient  mounting/interface  attachment  provisions. 

A  literature  survey  found  only  one  rotary  thermal  coupling  design  report'11,  the  abstract 
of  which  states: 

"This  report  documents  the  work  to  design,  fabricate  and  test  a  rotating 
thermal  joint  for  surveillance  satellite  applications.  The  rotating  thermal  joint 
uses  heat  pipe  technology  to  transfer  heat  into  and  out  of  the  joint.  Heat  is 
transferred  through  a  liquid  NaK  filled  annulus  between  the  rotating  and 
stationary  portions  of  the  joint.  The  joint  has  been  shown  to  transfer  1500 
watts  with  a  15C  temperature  drop." 

i 

This  rotary  thermal  joint  design  is  not  useful  in  the  cryogenic  temperature  range  due  to 
the  use  of  the  liquid  metal.  There  is  also  a  question  of  service  life  and  reliability  due  to 
the  requirement  for  a  rotating  mechanical  seal  to  retain  the  liquid  metal. 

The  design,  fabrication,  testing,  installation,  and  unsuccessful  launch  of  a  mission 
specific  rotary  thermal  coupler  has  been  reported121,  but  documentation  is  unavailable. 

12  Basic  Concept 

Consideration  of  the  general  requirements  for  a  rotary  thermal  coupling  listed  above 

suggests  that  the  only  truly  compatible  design  must  be  based  on  radiant  heat  transfer 

between  the  rotating  and  stationary  members.  Given  this  design  constraint,  it  is 
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apparent  that  there  are  only  two  possible  geometrical  variations;  concentric  cylinders 
and  parallel  disks.  Figure  1  depicts  these  variations. 


Figure  1  -  Rotary  Thermal  Coupling  Alternative  Geometries 
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Figures  2  and  3  show  isometric  views  detailing  these  alternate  geometries 


Figure  2  -  Rotary  Disk  Thermal  Coupling  Isometric 


Heat  Source  Side  Heat  Sink  Side 

-30  Degree  View  +30  Degree  View 


Figure  3  -  Rotary  Concentric  Cylinder  Thermal  Coupling  Isometric 
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To  implement  these  coupling  design  concepts  into  a  practical  device,  it  is  necessary  to 
incorporate  bearings,  insulation,  structural  supports,  housing,  mounting,  and  interface 
provisions.  For  most  installations  the  mounting  interface  will  be  a  Bearing  and  Power 
Transfer  Assembly  (BAPTA).  A  BAPTA  provides  the  basic  structural,  rotary  support, 
rotary  drive,  power  transfer,  and  control  and  data  transfer  between  stationary  and  rotary 
segments  of  a  space  vehicle.  BAPTA  technology  is  well  developed,  with  sources  for 
providing  “off-the-shelf  “state-of-the-art’’  hardware.  The  basic  BAPTA  configuration 
consists  of  the  following  elements: 

outer  stationary  structural  housing: 

i 

provisions  for  mounting  the  stationary  segment  of  the  space  vehicle, 
outer  bearing  race  housings, 
drive  motor  stator  mounting, 
power  take-off  stationary  segment  mounting, 
control  and  data  assembly  stationary  segment  mounting; 
inner  rotary  structural  housing: 

provisions  for  mounting  the  rotating  segment  of  the  space  vehicle, 
inner  bearing  race  housings, 
drive  motor  rotor  mounting, 
power  take-off  rotary  segment  mounting, 
control  and  data  assembly  rotary  segment  mounting; 
large  diameter  precision  bearings; 
large  diameter,  hollow  shaft  permanent  magnet  motor; 
rotary  transformer  or  slip  ring  power  transfer  subassembly; 
control  and  data  transfer  assembly  (optica!  or  electromagnetic); 
electrical  connectors; 

Figure  4  presents  a  general  schematic  of  a  BAPTA. 
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It  is  apparent  that  the  BAPTA  bearing  support  could  be  used  also  for  the  rotary  thermal 
coupling  bearing  support  in  a  highly  integrated  design.  It  is  also  apparent  that  provisions 
for  compensating  for  alignment  tolerances  will  be  required  for  designs  which  incorporate 
the  rotary  thermal  coupler  as  an  appendaged  subassembly  to  the  BAPTA.  The  details 
of  the  bearing,  support  and  thermal  insulation  design  issues  associated  with  a  highly 
integrated  BAPTA/Rotary  Thermal  Coupler  are  not  addressed  in  this  Phase  I  study  nor 
are  the  alignment  tolerance  compensation  issues  associated  with  the  appendaged 
subassembly  approach.  The  Phase  I  study  concentrates  on  the  design  issues 
associated  with  the  basic  rotary  thermal  coupler  heat  transfer  module. 

For  the  appendaged  subassembly  design  approach,  the  issue  of  bearing  integration 
with  the  basic  rotary  thermal  coupler  module  warrants  some  preliminary  consideration  to 
establish  the  potential  impact  of  these  bearing  systems  on  the  basic  rotary  thermal 
coupler  heat  transfer  module.  There  are  two  variations,  a  cold  bearing  design,  and  an 
ambient  bearing  design.  Figures  5  and  6  show  representative  implementations  using 
"cold"  bearings  while  Figures  7  and  8  show  representative  implementations  using 
ambient  temperature  bdarings.  It  is  evident  from  these  figures  that  the  rotary  disk 
configuration  has  intrinsic  design  advantages  over  the  rotary  concentric  cylinder 
configuration.  The  only  case  where  the  rotary  concentric  cylinder  approach  would  be 
the  configuration  of  choice  would  be  if  there  were  an  advantage  to  a  variable  area 
coupling  feature.  By  displacing  the  one  set  of  cylinders  axially  from  the  second  set,  the 
coupling  radiative  area  would  change.  It  is  possible  to  integrate  a  linear  positioning 
device  into  the  design  so  that  the  radiative  coupling  of  the  rotary  concentric  cylinder 
design  could  be  controlled,  thus  providing  a  method  of  temperature  control.  This 
approach  will  not  be  explored  as  part  of  this  research  and  development  effort  as  it  is 
beyond  the  scope  of  this  SBIR  program 
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Figure  6  -  “Cold  Bearing”  Rotary  Concentric  Cylinder  Design 
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Figure  7  -  Ambient  Bearing  Rotary  Disk  Design 


Figure  8  -  Ambient  Bearing  Rotary  Concentric  Cylinder  Design 
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The  "cold"  bearing  design  exhibits  some  mechanical  assembly  advantages  over  the 
ambient  bearing  design,  however,  the  "cold"  bearing  will  necessarily  be  of  the  dry 
lubricated  type  (lead,  molybdenum  disulfide,  etc.)13,4,51  while  the  ambient  bearing  design 
can  be  of  the  oil  lubricated  type  (Barden  precision  bearing,  ABEC  7,  Bray  oil  lubricated, 
shielded  metal  retainer).  The  "cold"  bearing  has  the  further  disadvantage  that  its 
frictional  energy  is  dissipated  within  the  cryogenic  coupling  path,  thus  adding  a  direct 
parasitic  heat  load.  It  does,  however,  introduce  an  additional  conductive  path.  On  the 
other  hand,  the  ambient  bearing  design  has  intrinsically  higher  structural  insulation 
parasitic  effects.  The  ultimate  choice  between  these  two  alternatives  requires  a  detailed 
design  analysis  which  will  be  undertaken  under  the  Phase  II  activities. 

Other  design  features  shown  in  Figures  5  through  8  are  the  polished  gold  plated 
radiation  shields,  .the  provisions  for  interfacing  with  the  rotating  heat  load  and  the 
stationary  heat  sink,  the  insulating  structural  mounting  provisions,  and  the  housing.  To 
minimize  the  parasitic  radiative  heat  load  from  the  ambient  housing  to  the  coupling 
elements,  the  interior  of  the  housing  should  be  polished  metal,  preferably  gold  plated, 
as  should  be  the  external  surface  of  the  stationary  element.  To  further  reduce  the 
radiative  heat  load,  a  single  radiative  shield  is  also  included  (at  least  in  the  initial  design 
iteration).  Consideration  of  the  radiation  equation  applicable  to  such  a  system,  which 
can  be  approximated  adequately  as  per  equation  1.0,  it  is  evident  that  the  addition  of  a 
single  radiative  shield  will  half  the  radiative  heat  leak. 
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The  rotating  heat  load  interface  is  designed  to  be  concentric  with  the  axis  of  rotation. 
The  stationary  heat  sink  interface  is  located  off  axis.  This  configuration  is  convenient 
from  a  design'  viewpoint,  and  further,  allows  for  the  provision  of  a  hole  concentric  with 
the  axis  of  rotation  that  could  be  adapted  to  laser  signal  transmission  from  the  rotating 
assembly  to  the  fixed  portion  of  the  installation.  Details  for  actual  method  of  attachment 
to  these  thermal  interfaces  will  be  established  during  the  Phase  II  research  and 
development  program.  Consideration  must  be  given  to  the  effects  of  differential 
expansion  of  the  axial  dimensions  of  these  interfaces  during  operation. 

The  insulating  structurail  mounting  provisions  for  the  rotary  cryogenic  thermal  coupler 
module  must  be  carefully  considered. 

The  parasitic  heat  leak  from  the 
ambient  housing  to  the  cryogenic 
coupler  elements  needs  ' to  be  kept  to 
a  small  fraction  of  the  heat  transmitted 
from  the  rotating  cold  heat  source. 

Techniques  which  can  accomplish 
this  include  pyramid  stacked  ball 
compression  supports,  low 
conductance  tension  supports,  and 
passive  orbital  disconnect  devices,. 

Figure  9  shows  the  first  .technique. 


Figure  9  -  Pyramid  Stacked  Ball  Support 

The  insulating  structural  support  must  be  capable  of  providing  for  the  ground  shock, 
handling,  and  the  launch  environment  loads  and  vibration,  must  provide  for  the 
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differential  thermal  expansion  between  the  housing  mounting  interface  and  the 
cryogenic  coupling  elements  during  operation,  and  must  not  introduce  a  large  parasitic 
heat  leak.  Any  of  these  techniques  can  be  successfully  adapted  to  the  Rotary  Cryogenic 
Thermal  Coupler  design,  the  appropriate  selection  being  dependent  upon  the  specific 
application.  The  pyramid  stacked  ball  compression  support  is  most  useful  in  very  low 
temperature,  small  size  applications;  the  low  conductance  tension  support  is  primarily 
intended  for  moderate  size  applications  throughout  the  total  temperature  range;  and  the 
passive  orbital  disconnect  devices  being  intended  for  heavy,  large  applications  at  the 
higher  cryogenic  temperature  range.  The  design  and  integration  of  the  detailed 
insulating  structural  supports  will  conducted  as  a  major  element  of  the  Phase  II  segment 
of  is  program. 

The  housing  provides  the  support  points  for  the  insulating  cryogenic  structural  support, 
the  mounting  interface  of  the  entire  unit  to  the  ambient  vehicle  structure,  bearing 
support  for  the  ambient  temperature  bearing  design  version,  possibly  provisions  for 
compensating  for  the  differential  thermal  expansion,  and  the  envelope  which  ties  all  the 
parts  together. 

The  rotary  cryogenic  thermal  coupling  elements  can  be  fabricated  of  copper,  aluminum, 
or  a  carbon  fiber  composite  -  high  conduction  material  being  mandatory  for  thermal 
efficiency  and  low  weight.  For  the  metallic  designs,  the  elements  can  be  fabricated  from 
bar  stock  using  plunge  tool  electrical  discharge  machining  to  cut  the  voids  which 
surround  the  fins,  or  can  be  fabricated  from  copper  pieces  by  diffusion  bonding  or  from 
aluminum  pieces  by  vacuum  brazing.  Soldering  is  not  considered  to  be  an  acceptable 
process  for  this  application  due  to  excessive  thermal  impedance  and  the  possibility  of 
unacceptable  thermally  induced  internal  stresses. 
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The  use  of  carbon  fiber  composites  would  require  the  application  of  adhesive  bonding 
techniques.  There  are  no  apparent  advantages  to  using  this  currently  labor  intensive, 
costly  process. 

To  keep  the  size  and  weight  low,  the  fins  and  the  gaps  should  be  kept  as  small  as 
possible.  The  diffusion  bonding  and  the  vacuum  brazing  techniques  require  larger 
tolerances  on  the  gaps  to  avoid  contact  due  to  the  possible  warping  of  the  parts  during 
the  high  temperature  processes. 

The  fin  thickness  is  constrained  by  both  manufacturing  limits  and  conduction  fin 
effectiveness.  It  will  be  necessary  to  conduct  a  through  thermal  analysis  of  the 
integrated  design  to  select  the  optimum  fin  thickness  for  any  given  design  point.  For 
very  low  temperature  cryogenic  couplers,  the  manufacturing  limits  and  provisions  for 
differential  thermal  expansion  are  anticipated  to  be  the  governing  factors;  for  higher 
temperature  applications,  the  fin  effectiveness  will  most  probably  be  the  controlling 
feature. 

Due  to  the  simplicity  of  this  design  and  the  need  for  high  design  margins,  it  was  initially 
proposed  to  use  the  linearized  version  of  the  radiation  equation  to  account  for  the  fin 
radiation  boundary  conditions.  In  this  manner,  the  classical  equations  for  fin 
effectiveness  in  a  convective  environment  could  be  used  rather  than  the  more  complex, 
radiation  boundary  conditions 

The  linearized  form  of  the  radiation  equation  is: 

j  =  nAfTa-TO  (2.0) 
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The  radiative  heat  transfer  coefficient  can  be  expressed  as: 


hr=ar.(T23  +  T22T1+T2Tl2+Tl3)  (2 .1) 

The  fin  effectiveness  is  a  function  of  hr,  the  material  thermal  conductivity,  the  fin  radius 
ratio,  the  fin  thickness,  and  the  radial  dimension,  which  can  be  shown  to  be  of  the  form: 


nf  = 


(2.2) 


From  the  literature  survey,  however,  it  was  discovered  that  this  linearization  would  result 
in  optimistic  performance  predictions  rather  than  conservative  predictions.  Chambers 
and  Sommers161  have  shown  that  the  use  of  the  linear  approximation  overstated  the  fin 
effectiveness  of  circular  fin  radiating  to  a  0  K  radiation  environment  by  as  much  as  60%. 
Thus,  an  analysis  of  the  rotary  coupling  module,  including  the  effect  of  the  4th  power 
radiative  heat  transfer  between  the  stationary  and  the  rotating  members  was 
undertaken.  A  complete  description  of  this  analysis  is  presented  in  Section  3. 


2.0  Phase  I  Objectives 


Table  1-1  through  1-3  summarize  the  Phase  I  effort  on  a  task  by  task  format. 
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3.0  Rotary  Thermal  Coupler  Module  Thermal  Analysis 

i 

The  governing  differential  equations  are  developed  from  fundamental  principles.  Figure 
10  is  an  isometric  view  of  a  set  of  rotary  thermal  coupler  module  elements. 


Figure  10  -  Module  Element  Set  Isometric 
Figure  11  shows  an  idealized  cross  section  through  such  a  typical  set  of  elements. 


Figure  11-  Module  Element  Set  Thermal  Characterization 
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The  heat  balance  on  the  elemental  differential  volumes  of  the  typical  pair  of  disks  can 
be  graphically  depicted  as  shown  in  Figure  12. 


Figure  12  -  Elemental  Differential  Volume  Heat  Balance  Depiction 


Based  on  the  model  shown  above,  the  basic  heat  balance  relationships  are: 

dTA 

qA=-2OT6k—^  (3.0) 


qB=-2w8k-^  (3.1) 

I 

^ab  -  _c1ba  =  oe^AB  2w(TA  -  Tg  )dr  (3.2) 


dqA 

dr 


dr  =  -2rc5k 


dr  =  qBA 


(3.3) 


•dr  =  — 27t5k  r 


tb  dTB 

-2  +  Hr  dr_CfAB 


The  analysis  assumes  that  TA  >  Te .  that  the  A  disk  is  attached  to  the  central  rotating 

shaft,  that  the  B  disk  is  attached  to  the  exterior  cooled  wall,  and  that  the  outer  diameter 

of  the  A  disk  and  the  inner  diameter  of  the  B  disk  are  effectively  insulated.  Note  also 
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that  8  is  V2  the  individual  disk  thickness,  since  the  central  plain  through  a  disk  is 
effectively  adiabatic. 


Consideration  of  the  thermal  boundary  conditions  imposed  on  a  typical  set  of  rotary 
thermal  module  elements  implies  that  the  temperature  distribution  through  such  a  set 
will  be  as  shown  in  Figure  1 3, 


Figure  13  -  Rotary  Thermal  Module  Element  Set  Temperature  Profile 


From  the  foregoing,  the  differential  equations  governing  the  thermal  steady  state 
response  of  a  pair  of  thermal  elements  can  be  expressed  as: 

d2  Ta  1dTA  cteFab  (Ta  -  Tr  ) 

— -f  +  -— T - ~  =  0  (4.0) 

d  r2  r  dr  k5 
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df.Tg_  IdTe  ct£^ba(Ta  ~ Tg  ) 
d  r2  +  r  dr  +  k8 

subject  to  the  following  boundary  conditions: 


Ta  =  TAi  @  r  =  n 


dje 

dr 


0  @  r  =  Tj 


dr 


=  0 


r  =  r„ 


Tb  =  TBo  @  r  -  r0 


Equations  4.0  and  4.1  can  be  transformed  to  a  more  convenient  form  by  using  the 
following  substitutions: 


0A  = 


Ta 

Tai’ 


Op 


R  = 


ro  -r, 


P  = 


The  resulting  normalized  equations  are: 
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d2  0A  1  d0A 
d  R2  +  Ri  1  dR 


ge  FABT^(ro-r,)2(0:-0g) 
kd 


=  0 


(5.0) 


d2©B,  1  d©B  TeFBATAj(r0  - ri)2(0A  - ©g)  n  „ 

W+: - - k8 - =  0  <5'1) 

r\  H - 1 

(p-1) 

Equations  5.0  and  5.1  can  be  further  simplified  by  introducing  the  dimensionless 
parameter: 


,  CT£FABTAi(r0  -r|)2 
k8 

The  form  of  the  reduced  normalized  equations  are: 

1  d©A 


(5.2) 


d2  0A 


d  R: 


R  + 


1  dR 


-M©a-0b)  =  o 


(P-1) 


(6.0) 


d2  0F 


■  +  ■ 


1 


d0f 


dR2'  1  dR 
R  (p-i) 


■+M0l-eS)  =  o 


(6.1) 


The  transformed  boundary  conditions  for  equations  6.0  and  6,1  are: 

0A  =  1  @  R  =  0 
d0B 


dR 


d©A 

dR 


=  0  (5)  R  =  0 


=  0(5)  R=  1 


0B  - 


'Bo 


1  Ai 


R  =  1 
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To  aid  in  solving  the  coupled  2nd  order  nonlinear  differential  equations,  6.0  and  6.1, 
the  following  transformations  can  be  used: 


d0‘  -Q 
dR 

(7.0) 

3  ® 

n 

-5 

(7.1) 

This  allows  equations  6.0  and  6.1  to  be  expressed  as 

Tr+  “i  -Wi-eS)-0 

r+(p-i> 

t 

(7.2) 

d 

1  +Mei-eS)  =  o 
r+(p“-^) 

(7.3) 

with  the  corresponding  boundary  conditions  being: 


©A  =1@  R  =  0 
VF  =  0  @  R  =  0 

0B  =  v2-  @  R  =  1 

1  Ai 


Q  =  0  @  R  =  1 

This  set  of  four  coupled  1st  order  nonlinear  differential  equations,  subject  to  the 
boundary  conditions,  can  be  solved  using  numerical  methods. 
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Note  that  from  equations  6.0  and  6.1  the  following  relationship  can  be  derived: 

42, 


d20, 


1  d©A  dz0B  1  d0B 


dR2  R  +  _l_  dR  dR2  R+  1  dR 


(P-1) 


(P-1) 


From  which: 


d(0A+0B)^  1  d(0A+0B)  ^  ^ 


dR" 


R  + 


1 

(P-1) 


dR 


(8.0) 


(8.1) 


By  defining  C  =  ©a  +  ©b  equation  8.1  simplifies  to: 

1  dC 


d2C 

dR2 


R  +  T  ... 

(p-1) 

dC 


1  dR 


=  0 


(9.0) 


Now  introducing  the  variable  O  =  — ,  equation  9.0  transforms  to 


dO  O 

which  results  in  the  relationship: 

1  do  dR 

O 


R  + 


The  solution  to  equation  10.1  is: 


O  =  Or 


which  inverse  transforms  into: 


1 

(P-1) 


1 


[R(p  - 1)  + 1] 


dC  _  dC 
dR  dR 


1 


o  [R(p-i)  +  i] 


(10.0) 


(10.1) 


(11.0) 


(ill) 
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and  thence  into: 


dQA  d0e  1 

dR  +  dR  l  dR  +  dR  JQ  [R(p  - 1)  + 1] 


(11.2) 


The  solution  to  equation  11.1  is: 


^  dC  ln[R(p  —  1)  + 1] 
4  4,0  dR  0  (p-1) 


(12.0) 


which  inverse  transforms  into: 


(0a+0b)  =  (0a+0b)o+(^-  ^  )ln[R(P~1)  +  1]. 

A  8/  A  60  dR  .  dR  '  (p-1) 


(12.1) 


Evaluating  equation  11.2  at  R=1  and  applying  the  appropriate  boundary  conditions: 


d0B  1  d0A 


(13.0) 


'  dR  |t  p  dR  |0 

which  is  simply  the  conservation  of  energy  boundary  condition.  Evaluating  equation 
12.1  at  R=1  and  applying  the  appropriate  boundary  conditions: 

/c\  d©A  I  ln(p)  ..  .  n. 


(0 A  +0b)i  “  (1  +  0b)o  +  ' 


(14.0) 


\  r\  O  '  l  \  O  '  U  jn  /  \  / 

dR  |„  (p-1) 

which  proves  to  useful  in  developing  a  numerical  analysis  convergence  criteria. 


3.1  -  Runge-Kutta  Numerical  Analysis 


The  solution  to  the  set  of  1st  order  nonlinear  differential  equations  7.0  -  7.3  can  be 
obtained  by  creating  a  vector,  VR,  evaluated  at  R: 
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dV  dV 

then  substituting  R  andVR  into  the  vector  — ,  to  obtain 


dR 


dV 

dR 


d0 , 


dR 

d0R 


dR 
d Q 

dR 

dT 

dR 


R,vr 


(16) 


having  obtained  —I 
y  dR1 


,  VR+AR  is  then  computed  using  the  4th  order  Runge-Kutta 


method. 


The  process  starts  at  R=0  and  proceeds  to  R=1  in  steps  of  AR.  The  process  is 
mathematically  outlined  below: 


Kn  ARx 


dV 

dR 


dV 


dR 


R.Vr 


V«,  vr  +  2 


R  <-  R  +  ■ 


AR 


K2  <-  AR  x 

dV 

<— 

dV 

dR 

K,,R 

dR 

Vk2 

k2 

<- VR  +^r 

R  2 

K,  <—  AR  x 

dV 

<— 

dV 

dR 

K,.R 

dR 
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_  AR 
R*-R  +  — 

vK,  <-  vR  +k3 


K4  4-  AR  x 


dV 

dR 


<— 


K,.R 


dV 


dR 


R.V, 


*3 


Vr+aR  ^R  + 


(K1+(2xK2)  +  (2xK3)  +  K4) 


The  computation  runs  in  (n+1)  steps  from  R=0  to  R=1  by  AR  =  — .  After  each  loop, 

the  results  at  R=1  are  compared  with  the  boundary  conditions,  the  error  computed, 
new  input  variables  at  R=0  are  then  generated  using  an  algorithm  designed  to  force 

i 

convergence  and  the  process  repeated  until  the  error  is  reduced  to  an  acceptable 
value. 


Given  the  temperature  distribution,  and  the  exact  boundary  conditions,  the  fin 
effectiveness  can  be  computed  from: 


-2K5- 


dTfi 


dr 


*f 


r=rt 


°e(r0  ~rj  )(Ta, 


tb40I) 


(17) 


which  reduces  to  the  form: 


~2Qr=q 

Mp-i)0-©br.,) 


(18) 


A  listing  of  the  APL  computer  program  which  was  written  to  compute  the  results  of  the 
analysis  described  above  is  provided  in  Appendix  I.  Several  cases  have  been 
computed,  the  results  of  which  are  presented  in  Appendix  II.  Using  these  results  Figure 
14  plots  r|f  versus  X  for  various  values  of  p  and  a  specific  value  of  0Br  i  . 
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0 

.2  .4  .6 

.  8 

1 

rc  (D0-D,)U,T= 

^  2  V  k5 

Figure  14  -  Fin  Effectiveness;  rif  Versus  X,  p 


The  computer  program  generates  a  complete  Rotary  Thermal  Coupler  design  using  the 
methodology  described  in  the  Section  4.0  Rotary  Cryogenic  Coupler  Preliminary 
Design. 
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4.0  Rotary  Cryogenic  Thermal  Coupler  Preliminary  Design 

The  preliminary  design  is  based  on  an  the  concept  of  an  interleaved  stack  of  radiantly 
thermally  coupled  disks,  such  as  shown  in  Figure  2.  Figure  15  defines  the  basic 
parameters  which  can  be  used  to  characterize  such  a  module. 
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Based  on  the  parameters  shown  in  Figure  15,  the  gap  between  an  adjacent  rotor-stator 
set  of  disks  is  g5.  thus  the  distance  between  the  adiabatic  planes  of  an  adjacent  rotor- 
stator  set  of  disks  is: 


s  =  (g  +  2)5 

Given  H,  the  overall  height  of  stacked  set  of  disks,  the  number  of  rotor-stator  radiant 
interfaces,  N,  in  a  module  is: 


s 


The  nominal  outside  diameter  of  the  thermally  radiating  area  of  a  disk  is  denoted  as  D0, 
the  nominal  inside  diameter  of  the  thermal  radiating  area  of  a  disk  is  denoted  as  Dj.  The 
effective  coupled  radiation  area  per  module  is  thus: 


.  "(D’-D’) 


N 


n 


n 


(d;-d?)h 

4(g  +  2)5 
d; 


i- 


D20) 


D.H 


4(g  +  2)5 


The  overall  diameter  of  a  module  is  greater  than  D0  by  twice  the  wall  thickness  of  the 
stator,  the  wall  thickness  being  established  by  the  conditions  imposed  due  to  the 
combined  influence  of  the  need  to  provide  a  low  thermal  impedance  conductive  path 
coupling  the  module  end  plates  with  the  stationary  disk  stack  together  with  the  annular 
area  needed  to  provide  sufficient  surface  area  to  assure  that  a  sound,  high  conductance 
metallurgically  bonded  stack  of  disks  can  be  fabricated. 
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The  overall  outside  diameter  of  a  complete  Rotary  Cryogenic  Thermal  Coupler  will  be 
greater  than  the  overall  outside  diameter  of  a  bonded  module  by  twice  the  radial 
dimension  required  by  the  thermal  support/insulation  system.  The  height  of  a  complete 
coupler  will  be  the  combined  height  of  the  module,  the  end  plates,  and  the  thermal 
support/insulation  system.  To  minimize  the  parasitic  heat  leak  into  the  module,  it  is 
desirable  to  maximize  the  internal  volume  per  unit  surface  area.  A  cylinder  with  the 
height  equal  to  the  diameter  has  the  maximum  volume  per  unit  surface  area.  Since  the 
thickness  of  the  thermal  support/insulation  system  plus  the  conductive  walls  are 
approximately  the  same  dimension  axially  as  radially,  the  stack  height  of  a  module,  H, 
excluding  end  plates,  should  be  approximately  equal  to  D0,  the  outside  diameter  of  the 
thermal  radiating  area  of  the  module  disks.  Thus  the  net  thermal  radiation  area  per 
module  is: 


7C 


V1'2' 


vp. 


D: 


Am  4(g  +  2)5 


(20.0) 


Based  on  the  fundamentals  of  radiation  heat  transfer  and  the  relationships  derived 
above,  the  governing  equation  for  the  heat  transfer/module  temperatures  is: 


Ql=n,Amasm(TL4-T4)  (21.0) 


In  conformance  with  the  rotary  thermal  coupling  module  geometry,  the  insulating 
housing  area  can  be  defined  as: 

7tD| 

As  =  kDsHs  +  2““  (22.0) 
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but  Ds  =  Hs  in  conformance  with  the  shape  of  the  rotary  thermal  coupling  module. 
Hence: 


Ac  = 


3tiD; 


(22.1) 


Recalling  equation  1.0,  the  governing  equation  for  module  housing  heat  leak  is: 


Qc  —  r 


oesAs(Ts4-TE4) 


1  + 


"AsN 

c  Aj:  > 


(1.0) 


(n  + 1) 


By  similar  reasoning  to  that  used  in  developing  equation  22,1,  it  is  apparent  that: 

ae  =  (23.0) 

The  identity  Te  =  Tm  is  also  evident,  thus  equation  1 .0  can  be  revised  to: 


Q< 


oesAs(Ts4-T^) 


1  + 


(d,Y 

<De, 


(1.1) 


(n  + 1) 


By  combining  equations  1.1,  20.0,  21.0,  22.1,  and  23.0,  the  ratio  of  the  heat  leak  to  the 
heat  load  can  be  expressed  as: 


Q< 


6ss(g  +  2)5D2s(Ts4-T^) 


Ql 


rifS 


m 


Dr 


1  + 


2^ 


(n  +  1) 


( 


(24.0) 


1- 


ypJ 


d03(t4-tl4) 


Electro  Thermo  Associates 
32 


Phase  I  Final  Report 

Rotary  Cryogenic  Thermal  Coupling 


By  denoting  ATU  =  Tm  -  TL  where  ATL  is  the  maximum  temperature  difference  across 
the  module,  equation  24.0  can  be  expressed  as: 


Q. 


6es(9  +  2)p.~  Ds 


f  (  T  \4 

•s 

~^"m  - 


-1 


1 

_J 

o 

f 

f  Ds^| 

2Y 

r 

|2> 

f 

r  AT^l 

4\ 

hfSm| 

1  + 

_ =L  i 

(n  +  1) 

1- 

D? 

1- 

V. 

vDEy 

J 

V 

J 

V 

^  '  m' 

J 

(24.1) 


Equation  24.1  is  of  interest  because  it  is  possible  to  introduce  some  order  of  magnitude 
approximations  and  obtain  a  quantitative  expression  for  the  relative  heat  leak.  The  order 
of  magnitude  approximations  are: 


—  Sr  1.25 

Dp 


De 

D, 


1.10 


-s-  *  1.38 

D0 

P  *  4 

^4*0.20 

Tm 

—  *  0.0025 

D0 


Electro  Thermo  Associates 
33 


Phase  I  Final  Report 
Rotary  Cryogenic  Thermal  Coupling 


(n  +  1)  <4 
(g  +  2)>  3 
rjf  =  0.85 


Thus  the  quantitative  expression  for  the  relative  heat  leak  is: 


Qs 

Qm 


>  0.000618 

v 


f  i  ^4 
‘s 


\ 

-1 


(25.0) 


For  a  typical  application,  Ts  =  250  K,  and  Tm  will  be  a  value  such  as  10  K,  35  K,  or  60  K, 
depending  on  the  specific  application.  The  corresponding  relative  heat  leaks  are  then: 


Q< 


Qr 

Q< 


Q 

Q< 


m 


Q, 


250K 

10K 

250K 

35K 

250K 

60K 


=  241 


=  1.61 

=  0.19 


It  is  apparent  that  the  relative  heat  leak  is  not  excessive  for  the  60  K  application,  and 
probably  not  excessive  for  the  35  K  application,  however,  for  the  10  K  application,  the 
relative  heat  leak  would  prove  to  be  unacceptable. 
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To  correct  this  situation,  an  active  thermal  shield  could  be  incorporated  within  the 
insulation  system.  The  active  heat  shield  could  be  heat  sunk  to  a  120  K  thermal 
radiator.  This  would  reduce  the  relative  heat  leak  for  these  application  to: 


250K 
Qs  T20K 

10K 

250K 

Qs  120K 

Qm 
m  35  K 

250K 

Qs  120K 

Qm 

m  60  K 


12.8 

0.085 

0.009 


The  reduction  in  relative  heat  leak  for  the  60  K  application  would  not  appear  to  be 
sufficient  to  justify  the  added  complexity  of  an  intermediate  120  K  active  heat  shield  ,  but 
it  could  be  justified  for  the  35  K  application.  The  relative  heat  leak  for  the  10  K 
application  would  probably  still  be  unacceptable. 

The  difficulty  with  the  10  K  application  could  be  resolved  by  incorporating  a  60  K  active 
heat  shield,  the  relative  heat  leak  would  reduce  to: 


Qjl 

Qm 


60  K 

10K 


0.80 


The  60  K  heat  sink  could  be  provided  by  the  first  stage  of  the  cryogenic  refrigerator 
required  to  provide  the  10  K  heat  sink  that  is  to  be  transferred  across  the  rotary 

1 

cryogenic  thermal  coupling. 
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To  reduce  the  heat  load  on  the  1st  stage  of  the  cryogenic  refrigerator,  the  addition  of  a 
120  K  active  thermal  shield  within  the  rotary  cryogenic  thermal  coupler  insulation 
system  should  be  considered. 


The  preceding  analysis  demonstrates  the  viability  of  the  concept  of  the  rotary  cryogenic 
thermal  coupler  from  a  thermal  management  viewpoint.  The  importance  of  an  adequate 
insulation  system  has  been  noted.  Fortunately,  the  technology  for  providing  such 
cryogenic  insulation  systems  is  state-of-the-art,  including  the  incorporation  of  actively 
cooled  shields  integrated  into  the  insulation  system.  To  provide  a  more  definitive 
mathematical  description  of  a  Rotary  Thermal  Coupler,  the  following  analysis  proves 
useful.  First,  rewrite  equation  5.2  as: 


X  = 


oeFABTL3D?(1---)2 

r 

4kS 


(5.3) 


Now  by  combining  equations  5.3,  20.0,  and  21.0,  assuming  that  FAB  =  1,  and  solving  for 
D0  one  obtains: 


(9  +  2) 


D„  = 


1-- 

v  py 


Ql 


r|f7rkX 
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ffl 
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f 

(j  } 

1- 

- 

1- 

m 
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^py 

y 

V 

ItlJ 

y 

(26.0) 


T, 


From  equation  26.0  and  the  definition  of  p ,  the  expression  for  Dj  is: 


D0 

D,  =  —  (27.0) 

'  P 

Having  computed  D0,  given  QL,  Tm,  TL,  p,  7,  nf,  a,  e,  g,  and  k,  it  is  possible  to 
compute  the  corresponding  value  for  5  from  equation  5.3  as: 
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f 


os\3D2 


1 


\2 


5  = 


1- 

V  PJ 


4kX 


(27.0) 


Based  on  the  geometry  presented  in  Figure  15  and  the  foregoing  definitions  and 
constraints,  the  number  of  rotor  disks  required  is: 


Do 

25(g  +  2) 


(28.0) 


and  the  number  of  stator  disks  is: 

Ns  =  Nr  +1  (28.1) 

Equations  20.0  and  21.0  can  be  used  to  verify  the  calculations  by  computing  Am,  and 
Ql.  The  relative  heat  leak  can  then  be  computed  by  using  equation  24.0.  Appendix  I 
presents  the  APL  computer  program  that  parametrically  performs  these  calculations. 
Appendix  II  presents  the  results  of  several  typical  analysis. 


In  conducting  these  calculations,  it  is  good  practice  to  use  a  value  for  QL  than  is  at  least 
50%  greater  than  the  actual  cooling  load  anticipated.  For  a  final  design,  It  is  necessary 
correct  the  value  of  5  to  correspond  to  14  available  stock  thickness.  The  actual  outside 
diameter  of  the  rotary  disk  elements  should  be  about  DB  =  D„  -45 ,  and  the  actual 

inside  diameter  of  the  stator  disks  should  be  about  DA  =  Dj  +45 .  The  inside  diameter 

of  the  rotor  disks  will  be  reduced  from  Dj  by  twice  the  wall  thickness  of  the  rotor 
conducting  shaft,  and  the  outside  diameter  of  the  stator  disks  will  be  increased  from  D0 
by  twice  the  wall  thickness  of  the  conducting  outside  structure  of  the  module. 
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5.0  Rotary  Thermal  Coupler  Module  Critical  Component  Demonstration 

It  has  been  concluded  that  the  critical  technology  that  needs  to  be  demonstrated  is  the 

i 

ability  to  manufacture  a  set  of  module  components  such  as  shown  in  Figure  2.  An 
obvious  method  to  be  considered  would  be  to  cut  the  stator  half  cylinders  from  bar  stock 
using  plunge  EDM  machining  and  the  rotor  using  wire  EDM  machining  employing  a 
moving  workpiece  adapter.  For  optimum  (minimum  weight)  design,  the  fins  should  be 
tapered.  Due  to  the  effect  of  radial  direction  of  heat  flow,  the  stationary  fins  attached  to 
the  interior  cylindrical  wall  can  be  thinner  than  the  rotating  fins  attached  to  the  central 
shaft  for  an  optimized  design. 

In  principle,  the  EDM  machining  approach  could  provide  the  capability  to  make  tapered 
fins.  Unfortunately,  the  tooling  for  the  EDM  process  would  be  very  expensive,  and  may 
not  be  able  to  provide  the  close  spacing  and  the  low  fin  thickness  to  fin  diameter  ratio 
that  is  required  to  meet  the  low  relative  heat  leak  goal. 

A  stacked,  bonded  sheet  metal  assembly  is  a  low  cost  method  for  demonstrating  the 
technology.  This  approach  can  provide  both  close  spacing  and  the  low  fin  thickness  to 
fin  diameter  ratio.  The  core  sizing  routine  outlined  in  Section  4.0  anticipated  this 
conclusion,  hence  the  reference  to  stock  material  thickness.  Several  alternative 
fabrication  techniques  were  evaluated.  The  basic  materials  considered  were  limited  to 
Aluminum  and  Copper.  Fabrication  techniques  considered  included  soldering,  brazing, 
diffusion  bonding,  and  spin  welding. 
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This  phase  of  the  task  consumed  many  man  hours,  requiring  an  extensive  review  of  the 
literature  and  numerous  conversations  with  potential  material  suppliers,  sheet  metal 
fabricators,  and  metal  bonding  job  shops. 


After  all  the  material,  processes,  and  vendors  were  surveyed,  it  was  concluded  that 
fabrication  of  the  rotary  thermal  module  elements  as  vacuum  brazed  Aluminum  parts 
was  the  best  choice.  For  this  demonstration,  a  small  module  was  designed  which 
incorporated  some  self  fixturing  provisions  to  simplify  the  fabrication  process.  Figures 
16  and  17  show  graphically  summarize  the  parts  details  and  their  assembly  process. 
The  parts  list  for  a  module  was: 


Rotary  Thermal  Coupler  Parts  List 


ITEM 

SIZE 

MATERIAL 

QUANTITY 

Rotor  Disk 

0.50”1D,  3,00"OD 

0.010”  3003  Al 

40 

Stator  Disk 

1.00’ID,  3.50"OD 

0,010”  3003  Al 

39 

Stator  End  Disks 

1.00”ID,  3.50”OD 

0.025”  3003  Al 

2 

Rotor  Spacer 

0.50"ID,  0.75"OD 

0.020"  #8  Braze  Sheet 

41 

Stator  Spacer 

3.10"OD,  3.50"OD 

0.020"  #8  Braze  Sheet 

40 

Rotor  Stack  Tube 

1.75”  Long 

0.50”OD,  0.10”  Wall  3003  Al 

1 

i 

The  parts  were  cut  from  the  sheet  stock  and  the  cut  edges  were  deburred.  The  rotor 
disks  and  spacers  were  assembled  on  the  rotor  stack  tube.  The  stator  end  disks,  stator 
spacers,  and  stator  disks  were  assembled  in  a  3  pin  assembly  braze  fixture,  the  pin 
spacing  defining  the  stator  outside  diameter. 
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Figure  16  -  Rotor  Parts  and  Assembly 
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Figure  17  -  Stator  Parts  and  Assembly 
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The  two  assemblies  were  placed  within  specially  designed  vacuum  brazing  fixtures, 
placed  within  the  vacuum  furnace,  thermocouples  were  installed  to  monitor  the 
temperature  response  of  the  parts  being  brazed,  and  then  the  parts  were  vacuum 
brazed  using  a  temperature-time  profile  proprietary  to  the  vacuum  brazing  vendor. 

After  the  successful  brazing,  the  stator  was  cut  diametrically  in  half.  The  two  stator  parts 
were  assembled  over  the  rotor  part.  Figure  18  presents  photographs  of  the  details  and 
the  fabricated  parts. 


Figure  18  -  Photograph  of  Details  and  Assemblies 
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6.0  Conclusions  and  Recommendations 

This  Phase  I  effort  has  developed  and  demonstrated  the  principles  and  practices 
required  for  the  design  and  fabrication  of  a  Rotary  Cryogenic  Thermal  Coupling. 

A  detailed  mathematical  analysis  of  the  radiant  heat  transfer  occurring  within  the  basic 
Rotary  Cryogenic  Thermal  Coupling  module  has  been  presented,  resulting  in  a 
mapping  of  the  module  disk  fin  effectiveness  versus  the  disk  radius  ratio  and  a 
component  dimensionless  characterization  parameter,  X ,  for  a  specific  temperature 
regime. 

An  APL  computer  program  has  been  written  to  parametriclly  size  a  Rotary  Thermal 
Coupler  based  on  the  mathematical  analysis  developed  in  this  Phase  I  effort.  A  listing  of 
the  program  is  presented  in  Appendix  I.  A  printout  of  the  results  of  several  different 
cases  are  presented  in  Appendix  II. 

The  thermal  analysis  has  demonstrated  that  the  thermal  heat  leak  parasitic  losses  can 
be  kept  within  a  reasonable  percentage  of  the  transferred  heat  by  incorporating  active 
thermal  shielding.  The  technology  is  applicable  down  to  10  K. 

The  preliminary  design  analysis  concluded  that  the  unproven,  critical  component  portion 
of  the  Rotary  Cryogenic  Thermal  Coupling  was  the  basic  radiation  coupled  disk 
assembly.  Several  alternate  materials  and  processes  approaches  were  considered.  It 
was  concluded  that  a  concept  based  on  vacuum  brazed  Aluminum  fabrication  was  the 
best  option. 
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A  proof  of  concept  module  was  designed  and  fabricated.  This  task  demonstrated  that 
the  module  components  could  be  fabricated.  For  the  purposes  of  obtaining  as  much 
information  about  this  design  approach,  the  fin  spacing  was  selected  to  be  at  the 
extreme  minimum  value,  namely  a  value  of  g  =  1  was  selected,  resulting  in  a  nominal 
clearance  of  0.005”.  This  value  is  now  considered  to  be  too  close  for  designs  using  3” 
diameter,  0.010”  thick  3003  Aluminum  since  a  slight  warpage  of  either  the  rotor  or  the 
stator  disks  will  result  in  interference.  While  the  brazed  parts  initially  assembled  with 
adequate  clearance  to  provide  for  frictionless  rotation,  subsequent  handling  has 
resulted  in  sufficient  distortion  of  the  soft  Aluminum  parts,  thus  causing  rotor-stator 
contact.  It  is  possible  that  this  distortion  could  be  removed  by  reforming  the  parts  using 
some  sort  of  tooling  to  straighten  out  the  disks.  It  is  recommended  that  at  this  time, 
however,  that  a  value  of  g  >  2  be  used  for  future  initial  designs. 

A  proposal  for  Phase  II  was  not  submitted  because  there  was  no  formal  request  issued 
by  the  contracting  agency.  The  features  which  should  be  considered  in  such  a  follow-on 
effort,  nevertheless,  were  considered.  A  partial  list  of  such  recommended  future  work 
should  include: 

1.  Development  of  an  integrated  computerized  design  program. 

A.  The  design  process  and  the  basic  kernel  of  the  fin  effectiveness 
program  could  be  melded  into  an  computer  program  that  would  output  a 
detailed  parts  list,  given  the  design  heat  load  and  temperature.  The 
program  should  include  an  internal  table  listing  the  available  Aluminum 
sheet  stock. 
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2.  Sources  and  stocks  of  Aluminum  vacuum  brazing  sheet  should  be  identified. 
A.  Although  Reynolds  Metals  Co.,  Richmond,  Virginia  manufactures 
Aluminum  vacuum  brazing  sheet,  primarily  for  the  automotive  industry, 
the  minimum  buy  requirement  is  so  great  as  to  virtually  eliminate  them  as 
a  potential  source.  An  intensive  search  failed  to  locate  any  stocking 
distributors.  The  material  used  in  the  fabrication  of  the  experimental  unit 
was  finally  obtained  as  an  addendum  to  an  order  from  a  Scandinavian 
supplier  to  the  vacuum  furnace  brazing  vendor.  This  is  not  a  satisfactory 
solution  for  flight  type  production  hardware. 


3.  A  pair  of  engineering  prototype  units  should  be  designed,  fabricated,  and 
tested.  Figure  19'  outlines  the  recommended  testing  procedure. 
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6.1  Potential  Technology  Use 

The  Rotary  Cryogenic  Thermal  Coupler  could  prove  to  be  a  cost  effective 
component  for  use  on  low  cost,  light  weight  surveillance  satellites  to  aid  in  the 
thermal  management  of  infrared  detector  systems,  optics,  and  superconducting 
electronics.  The  technology  is  adaptable  to  higher  temperature  applications,  such 
as  vehicle  electronics  temperature  control  and  deployable  radiator  systems  when  it 
is  desirable  to  transfer  thermal  energy  across  a  mechanical  joint.  A  pair  of  couplers 
can  be  assembled  to  provide  effectively  a  "universal  joint"  thermal  path,  thus  raising 
the  possibility  of  eliminating  the  need  to  biannually  "flip"  vehicles  which  incorporate 
deep  space  viewing  cryogenic  thermal  radiators. 

Rotary  Cryogenic  Thermal  Couplers  could  be  employed  in  the  medical  industry  in 
MRI  scanning  systems  and  cryogenic  surgery,  Such  a  device  could  be  used  to 
enhance  the  application  of  superconducting  devices  by  allowing  the  relative  rotary 
movement  between  the  cryogenic  heat  sink  (refrigerator  or  cryogen)  and  the 
cryogenic  heat  source  (superconducting  device).  It  could  be  employed  in  SQUIB  low 
intensity  magnetic  field  scanning  equipment. 


Electro  Thermo  Associates 
46 


Phase  I  Final  Report 

Rotary  Cryogenic  Thermal  Coupling 


7.0  References 

1.  Shaubach,  Robert  M.;  "Rotating  Heat  Pipe  Joint  -  Design,  Fabrication,  and  Test"; 
Thermacore,  Inc.;  AFWAL-TR-83-3053;  June  1983 

2.  Backovsky,  Z.F.,  Rockwell  International,  Personal  Communication, 

3.  McConnell,  B.D.,  Mecklenburg,  K.R.;  "Unlubricated  Bearing  Operation"; 
Conference  on  Wear  and  Lubrication;  London;  September  1967 

4.  Mecklenburg,  K.R.;  "Performance  of  Ball  Bearings  in  Air  and  Vacuum  With  No 
Added  Lubrication";  Midwest  Research  Institute;  Air  Force  Contract  No.  F33615- 
69-C-1236;  December  1971 

5.  Todd,  M.J.;  "Solid  Lubrication  of  Ball  Bearings  for  Spacecraft  Mechanisms"; 
Tribology  International,  December  1982 

6.  Lieblein,  S.;  "Analysis  of  Temperature  Distribution  and  Radiant  Heat  Transfer 
Along  a  Rectangular  Fin  of  Constant  Thickness";  NASA  Technical  Note  D-196; 
1959 


Electro  Thermo  Associates 
47 


Phase  I  Final  Report 
Rotary  Cryogenic  Thermal  Coupling 


8.0  Appendix  I 

[ 0 ]  RTC 

[1]  a  Program  sizes  a  Rotary  Thermal  Coupler  Module 

[2]  a  The  process  is  based  on  that  outlined  in  the 

[3]  a  SBIR  Phase  I  Final  Report  "  Rotary  Cryogenic 

[4]  a  Thermal  Coupler",  Electro  Thermo  Associates 

[5]  a  31  October  1896,  Dr.  A.L.  Johnson,  P.I. 

[61  'RTCInput 1 

t7]  RTCInput 

[8]  'RTCInitalization ' 

[9]  RTCInitalization 

[10]  'RTCAnalysis ' 

[11]  RTCAnalysis 

[12]  1 RTCOutput 1 

[13]  RTCOutput 


[0]  RTCInput 

[1]  a  Function  is  the  input  module  for  the  RTC  module  sizing  program.  Note  that 

[2]  a  Copyright,  Dr. A.L.  Johnson,  1996 

[3]  'Input  module  for  the  Rotary  Thermal  Coupler  design  program.  Note  that  the1 

[4]  'program  accepts  vector  strings  for  the  Radius  Ratio  and  the  Lamda  input' 

[5]  'terms,  which  expedits  parametric  studies.  The  terms  in  the  vector  strings' 

[6]  'are  to  be  seperated  by  spacebar  blanks.  Due  to  printout  limitations,  the  vector' 

[7]  'length  of  the  Radius  Ratio  parameter  is  limited  to  6  quantities.' 

[6] 

1 9  ]  L1:QL*-aIN  'Enter  the  design  heat  load,  including  margin,  0.01  <  QL  <  20  watts; 
[10]  -1 <QL<0.01)v(QL>20) ) /Li¬ 

lli]  L2 : TL-aIN  'Enter  the  heat  source  temperature,  at  the  RTC  junction,  5K  <  TL  <  500K 

[12]  -( (TL<5 )v(TL>500) )/L2 

[13]  L3:Tm*-JtIN  'Enter  the  heat  sink  temperature  at  the  RTC  interface,  4K  <  Tm  <  (Tl-1) 

[14]  -{ <Tm<4 )v(Tm>(TL-l) ) )/L3 

[15]  LA: Tamb~.iIN  'Enter  the  ambient  environment  temperature,  TL  <  Tamb  <  BOOK 

[16]  -( (Tamb<TL)v(Tamb>500) )/LA 

[17]  a  L4:a‘-aIN  'Enter  the  disk  thickness,  a,  in  mm,  0.1mm  <  a  <  lmm;  ' 

[18]  a  -(  U<0.1)vU>l)  )/L4 

[19]  L5:K*-xIN  'Enter  the  disk  thermal  conductivity,  k  -  watts/cm. K;  ' 

[20]  -( (K<0.001)v(K>100) )/L5 

[21]  L6:g*-AIN  'Enter  the  disk  gap  factor,  g,  where  gA=gap,  1  <  g  <  5;  ' 

[22]  -( (g<l)v(g>5) )/L6 

[23]  L7  :mv*-  (  1 0 )  ,mv*-AlN  'Enter  the  Module  Radius  Ratio  Vector,  1.5  <  P  =  DO/Di  <  10; 

[24]  -(+/( (mv<l. 01) v(mv>15j ) )/L7 

[25]  -( (Pmv)>6 )/L0 

[26]  L8  :Lv-(  10 )  ,Lv-i.IN  'Enter  the  Lamda  parameter  vector,  0.0001  <  i  <  1  -  ' 

[27]  -(+/( ( Lv<  0 . 0  0  01  ) v ( Lv> 1 ) )  1/L8 

[28]  a  L9 :n-iIN  'Enter  the  number  of  numerical  integration  steps,  10  <  n  <  1000;  ' 

[29]  a  -( (n<io)v(n>1000) )/L9 

[30]  -End 

[31]  L0: 'Program  Units  the  allowable  number  of  terms  in  the  P  vector  to  6' 

[32]  -L7 

[33]  End : 
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[0]  RTCInitalization;X ; Y 

[1]  a  Function  establishes  the  initial  conditions  for  the  RTC  analysis 

[2]  n-300 

[3]  sigma*-5  .  67E“12  A'Steffan  Boltzman  constant,  Watts/sq  .  cm-K*4 

[4]  TA1-TL  A  K 

1 5  ]  m*-mv  [  1  ] 

[6]  em<-0.8  5 

[7]  esem«-0.03 

[8]  DsDo-1.25 

[  9  ]  TsTm*-Tamb+Tm 

[10]  DsDe-1.1, 

[11]  L-Lv[l] 

[12]  X-.01, ( 0.1x( 110) ) 

[13 ]  Y-0 .475 ,0.5125,0.5625,0.65,0.76,0.89,1.035,1.185,1.34,1.51,1.65 

[14]  Co^YSX • . * ( ( IPX ) -1 )  o  a  coefficients  for  i  power  series  for  UOref 

[15]  R-0 

[16]  Slref «-S0«-TnH-TL  o  a  establishes  initial  SO  and  final  SI  values 

[17]  T0-1 

[18]  U0*-U0ref 

[19]  VO-O 

[2  0]  QsQlv*-Div*-Dov*-Av-Amv<-QLv<-Nv<-W0*-WN*-Nf  «-M«- 1  0 

[21]  Wio-Wn-S0,T0,U0,V0 

[22]  dR-l+n 

[23]  k-1 

[24]  yi 


[0]  Result*-U0ref  ;  a  ;b ;  c 

[1]  SlTl-2xS0ref * ( +/Co* ( L*0 . 5 ) * (  (  ipCo )-l )  ) 

(2  ]  Result<-_10E~6  L  (m-l)x(siTl-(  1+S0ref  >  )  -f-  ( <*>m ) 


[0]  RungeKutta ;K1 ;K2 ;K3 ;K4 

[1]  a  Function  runs  4th  order  Runge-Kutta  numerical  analysis 

[ 2 ]  W-Wn 

[3]  Kl-dRxdWdR 

[4]  R-R+(dR*2) 

[5]  W-Wn+(K1^2) 

[6]  K2-dRxdWdR 

[7]  W*-Wn+  ( K2t2  ) 

[8]  K3-dRxdWdR 

[9]  R«-R+  ( dR-r2  ) 

[10]  W<-Wn+K3 

[11]  K4~dRxdWdR 

[  12  ]  Wn-Wn+(  (Kl+(  2xK2  )  +  (  2xK3  )  +K4  )  -=-  6  > 
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[ 0 ]  RTCAnalysls ; cf ; nf ; Do ; a a ; Nn ; Am ; q 1 ; QsQ 1 ; A; B 

[I]  fi  Function  solves  the  Rotary  Thermal  Coupler  differential  equations 

[21  '  l ,  lamda;  P*=Do/Di;  Iterations' 

[3]  L2:M-M,Wn 

[4]  RungeKutta  o  a  Function  runs  4th  order  RungeKutta  numerical  analysis 
CS]  -(R<l)/L2 

[6]  -((WnC4)<0)A(((|  (m-  ( UO+Wn  [  4  ]  )  )  )xiOE5  )  <5  )  )/L3 

[7]  U0-(U0-( (Wn[3]x2+(1?5) ) ) ) 

[8]  cf-(l?9)+20 

[9]  U0-(U0+(cfx<  (m-1  )T®m)  x  (  (Wn[l]+Wn[2]  )-(l+S0)  )  )  l-Ml+cf  ) 

[10]  R-0 

[II]  Wn-S0,T0,U0,V0 

[12]  Wi-Wn  0  M-tO 

[13]  j-j  +  1 

[14]  •»  (  j>8  00  )  /L7 

[15]  -»L2 

[16]  L3 : ASO-Slref-Wntl] 

[17]  aS0-(1-(L+(1+(1?5)  )  )  )xaS0 

[18]  -<  (  I  USO  +  Slref  )  )<10E~5  )/L4 

[19]  S0-S0+AS0 

[20]  •*(  (S0<0)v(S0>l)  )/L6 

[21]  -*L2 

[22]  L4 :Wi-M[l] ,M[2] ,M[3] ,M[4 ] 

[23]  M-M,Wn  0  W0-W0,W1 

[24]  WN-WN,Wn 

[25]  Nf-Nf  ,nf— (  (2xWi[3]  )*  (Lx  (m+1 )  x  ( l- (Wn[  1  ]  *4  )  )  )  ) 

[26]  Do-(QLx(g+2  )x(  (l-(l+m)  )*2)  )  +  (  nf  xTLx  ( l- (  (Tm+TL)*4  )  )xoix(i-(  (l+m)*2)  )xKxL) 

[27]  Dov-Dov,Do  0  Dlv-Div,  (Do-fm) 

[  28  ]  4v-av,44<-(  slgmaxemx  ( TL*3 ) x ( Do*2 ) x (  ( 1-  ( 1+m)  )  *2  )  )  +  ( 4  xRxL ) 

[29  ]  Amv-Amv,  Am-  (  olx  ( l  -  (  ( 1+m)  *2  )  )  x ( Do*3 )  )+(4x(g+2)xAA) 

[30]  QLv-QLv ,  q  1  -  ( nf  xAmxs  lgmaxemx  (  (  TL  *  4  )  -  (  Tm*  4  )  )  ) 

[31]  Nv-Nv,Nn-Do+  (  (g  +  2  )  xaa  ) 

[  32  ]  QsQl-6xesemx  (g  +  2 ) xaax  (  ( DsDo ) *2 ) x (  (  ( TsTm ) *  4  )  —  1  )  +Do 

[33]  QsQl-QsQl+ (nf x (  l+(DsDe*2  )  ) x4  x ( l- (  (  1-rm)  *2  )  )x  (  (  (  TL+Tm )  *4  )  - 1  )  ) 

[34]  QsQlv-QsQlv.QsQl 

[35]  (8  6TL) , (8  3Tm) , (8  OTj ) 

[36]  L6  :M-i0  0  j-1  o  i-LviL 

[37]  •*  (  i*=PLv )  /L7 

[38]  i-i  +  1 

[39]  L-Lv[  i  ] 

[40]  UO-UOref  o  SO-Slref 

[41]  Wn-SO ,T0 ,U0 , V0 

[42]  -»L2 

[43]  L7 : i-mvlm 

[4  4]  •»(  l=Pmv)/End 

[45]  i-i+1 

[46]  m-mv[l]  0  L-Lv[l] 

[47]  UO-UOref  P  SO-Slref 

[48]  Wn-SO , TO ,U0  ,  V0 

[49]  -L2 

[50]  End : 
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to] 

til 

[21 

[3] 

[41 

151 

[61 

[7] 

[8] 
[91 
tiol 
[111 


Result-dWdR;dSdR;dTdR;dUdR;dVdR;S;T;U;V  the 

o  The  der ivitlves  defined  above  are  used  to  compute 

s-orwtll  A  T-0rW[2]  O  U<-0LW[3]  o  V-0LW141 
dSdR*-V 

dUdR*-(Lx(  (T*4  >  —  (  S *  4  )  )  )  -  ( U-MR+ ( 1*  (m-1 )  )  )  ) 
dVdR<-ri)x(  (LX(  (T*4  )-(S*4  )  )  )  +  |VHR+(lr(m-l)  )  )  ) 

Resu  1 t -dSdR ,  dTdR ,  dUdR ,  d  VdR 


[0]  RTCOutput 

[1]  InputDataFile 

[2]  FinEff ectlvenessOutput 

[3]  OutsldeDlameterOutput 

[4]  InsldeDlameterOutput 

[5]  DlskThicknessOutput 

[61  HeatLeakParameterOutput 

[7]  RotorStatorDiskSpacerOutput 

[8]  RotorDlskCountOutput 
[S]  StatorDlskCountOutput 
[10]  ClearWorkSpaceVarlables 


[0]  InputDataFilejA 
[11  ’  ’ 

[2]  A«- 1  Heat  Load  ',(4  2TQL),'  watts;  Source  temp.,  TL ,  ’ 

[3]  A«-A,(6  2TTL),'  K;  Sink  temp.,  Tm,’,(5  2TTm),'K' 

[4  1  A 

[5] 

[6]  A*-' Ambient  Temp.f  Tamb,'(6  2TTamb) 

[7]  A«-A,  'K;  Thermal  conductivity,  k,’,(5  2  ▼K )  ,  '  watts/cm. K' 

[8]  A 

[9]  ’  ' 

[10]  A«-'  Disk  Gap  Factor,  g  =  ',(3  lrg) 

[11]  A*-A,  '  ;  Number  of  numerical  intergration  steps,  n  “',(5  OTn) 

[12]  A 

[13]  ’  ' 


[01 

[1] 

[2] 

[3] 


ClearWorkSpaceVarlables 

□ERASE  'Am  DsDe  Nf  QsQlv  TAi  Ustart  Wio  g 
□ERASE  'Wn  1  nf  av  Co  K  Nv  50  Tamb  W  cf  j 
□ERASE  'Lv  QLv  Slref  TsTm  WN  em  m  a  dR  k 


n  at  Amv  DsDo  Nn  R  TL  VO  esem' 
q 1  a  a  Div  L  QL  S1T1  Tm  WO  mv ' 
sigma  Do  Dov  M  QsQl  TO  U0  W1  as 
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[o] 
[i] 
12] 
C  3 1 

[4] 

[5] 
[61 
17] 
[8] 
[9] 


FlnEf f ect ivenessQutpu t ; Column ; Space ; Row ; Data ; output 


Column-10  3Tmv 

Space*-!'  i  \  '  )  -  (  (PColumn  )  P  1 
Column*- (  '  \  P  '  ) -Column 

Row*-(  (PLv)  ,  7  )  P  (  7  5  tLv  ) 

Data*-®  (  (  ( pmv )  ,  ( PLv )  )  PNf  ) 
output-RoWy [ 2 1 ( 10  3TData) 

'  p  *=  rO/ri,  lambda 


versus  radiation  fin  effectiveness' 


Columnr (Space-output ) 


[  0 ]  Ins ideD lame t erOu  tpu  t ; Column ; Spa  ce ; Row ; Da  t  a ; ou  tpu  t 

[ll  Column*-10  3Tmv 

[2]  Space*- (  '  i  \  ')  -  (( PColumn  ) P  ’  ') 

[3]  Column*-!'  \  p')yColumn 

[4]  Row*-!  (PLv),  7  )  p  (  7  5TLv) 

[51  Data*-®  ((  (Pmv)  ,  (PLv)  )P  (Dlv)  ) 

[6]  output-Rowy [ 2 1 ( 10  3TData) 

[7] 

[81  '  P  “'rO/rl,  lambda  -  i ,  versus  radiation  disk  ID  -  cm.' 

[91  '  ' 

[  1 0 1  Columny ( Spaceyoutput ) 


[0]  Ou ts IdeD iameterOu tpu t ; Column ; Space ; Row; Data ; output 

[1]  Column *-10  3Tmv 

[2]  Space*- (  '  i  \ ')  T  (( PColumn  )  P  ’  '  ) 

[3]  Column- ( 1  \  P')TColumn 

[4]  Row- ((PLv),  7  ) P ( 7  5 TLv) 

[51  Data-® (( (Pmv) , (PLv) )PDov) 

[6]  output-Rowr [2 1 ( 10  3TData) 

[71 

[81  '  P  =  rO/ri,  lambda  -  i,  versus  radiation  disk  OD  -  cm.' 

[91 

[101  Columny ( Spaceyoutput ) 


[0]  DlskThlcknessOutput , -Column ; Space ;Row;Data ; output 

[1]  Column-10  3Tmv 

[21  Space- ( ’  i  \ ’) - (( PColumn ) P ’  ') 

[31  Column-!'  \  P')y Column 
[4  1  Row- ((PLv),  7  ) P ( 7  5 TLv) 

[5  1  Data-® ( ( (Pmv) , ( PLv) )P ( 2  0*av) ) 

[61  output-Rowy [ 2 1 ( 1 0  3TData) 

[71  '  ’ 

[ 8 1  '  P  =  rO/ri,  lambda  -  i,  versus  radiation  disk  thickness  -  mm 

(91 

[10]  Columnr ( Spaceyoutput ) 
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[0]  HeatLeakParameterOutput ; Column ; Space ; Row ; Data ; output 

[1]  Column*-10  3Tmv 

[2]  Space*- (  '  i  \  '  ) 7  (  < PColumn ) P  1  '  ) 

[3]  Column*- (  '  \  P  '  ) yColumn 

[4]  Row*- (  ( PLv )  ,  7)P!7  5TLv) 

[5]  Data*-$  (  (  (Pmv)  ,  (P^v)  )P(QsQlv)  ) 

[6]  output*-Row7  [2  ]  ( 10  4TData) 

1  p  »  rO/ri,  lambda  —  x,  versus  Heat  Leak  Parameter,  Qs/QL 

[9] 

[10]  Column- ( Spaceyoutput ) 


[  0 ]  RotorStatorDiskSpacerOutput ; Column ; Space ; Row; Data ; output 

[1]  Column*-10  3Tmv 

[2]  Space*- (  '  x  \  1  )y  (  (PColumn)P  '  ’) 

[3] -  Column*- (  '  \  p')-Column 

[4]  Row*-(  (PLv)  ,  7  ) P  (  7  5TLv) 

[  5  ]  Data*-6)  (  (  ( pmv)  ,  (  PLv)  )P  (  2  Oxavx  ( 1+g  )  )  ) 

[6]  output*-Row- [2  ]  ( 10  3TData) 

[7]  '  ' 

[8]  '  p  =  rO/ri,  lambda  -  1,  versus  Rotor/Stator  Stack  Spacing 

[9]  '  • 

[10]  Columnr ( Spaceroutput ) 


[  0 ]  RotorDlskCountOutput ; Column ; Space ; Row; Data ; output 

[1]  Column*-10  3Tmv 

[2  ]  Space*- (  '  x  \ ') 7  (  (PColumn )P  '  ’) 

[3]  Column*- (  '  \  p;)TColumn 

[4]  Row*- ((PLv),  7  ) P  (  7  5 TLv ) 

[5]  Data*-b(  (  (pmv)  ,  (PLv)  )p(  r  (Nv-r2  )  )  ) 

[6]  output*-Row7  [2  ]  ( 10  OTData) 

[7]  1  ' 

[8]  p  =  rO/ri,  lambda  -  1,  versus  number  of  rotor  disks  ' 

[9]  ’  ' 

[10]  Column7  (  Spacerout'put ) 


[  0 ]  StatorDlskCountOutput ; Column ; Space ; Row; Data ; output 

[1]  Column*-10  3Tmv 

[2]  Space*- (  '  1  \  ')  7  (( PColumn  ) P  ’  ') 

[3]  Column*- (  '  \  P' ^Column 

[4]  Row*- (  ( PLv )  ,  1)9(1  5 TLv ) 

[5  ]'  Data*-0  (  (  (Pmv)  ,  (PLv)  )p  ( l  +  r  ( Nv-r2  )  )  ) 

[6]  output*-RowT  [  2  ]  ( 10  OTData) 

[7]  •  ’ 

[8]  p  =  rO/ri,  lambda  -  1,  versus  number  of  Stator  disks  ' 

[9] 

[10]  Columnr ( Space70utput ) 


mm ' 


Electro  Thermo  Associates 

53 


Phase  I  Final  Report 
Rotary  Cryogenic  Thermal  Coupling 


9.0  Appendix  II 

RTC 

RTC Input 

Input  module  for  the  Rotary  Thermal  Coupler  design  program.  Note  that  the 
program  accepts  vector  strings  for  the  Radius  Ratio  and  the  Lamda  input 
terms,  which  expedlts  parametric  studies.  The  terms  in  the  vector  strings 
are  to  be  seperated  by  commas.  Due  to  printout  limitations,  the  vector 
length  of  the  Radius  Ratio  parameter  is  limited  to  6  quantities. 


Enter  the  design  heat  load,  Including  margin,  0.01  <  QL  <  20  watts;  1 

Enter  the  heat  source  temperature,  at  the  RTC  junction,  5K  <  TL  <  500K;  65 

Enter  the  heat  sink  temperature  at  the  RTC  interface,  4K  <  Tm  <  (Tl-1);  GO 

Enter  the  ambient  environment  temperature,  TL  <  Tamb  <  500K  250 

Enter  the  disk  thermal  conductivity,  k  -  watts/cm. K;  2.5 

Enter  the  disk  gap  factor,  g,  where  gA=gap,  1  <  g  <  5;  2 

Enter  the  Module  Radius  Ratio  Vector,  1.5  <  P  =  D0/D1  <  10;  1.5,3,6,10 

Enter  the  Lamda  parameter  vector,  0.0001  <  I  <  1  -  ( . 01, . 02 , . 05 , . 1, . 15 , ( . l+( . Ixi9 ) ) )*2 


RTCInltallzation 
RTCAna lysis 

Heat  Load  1.00  watts;  Source  temp.,  TL,  65.00  K;  Sink  temp.,  Tm,60.00K 
Ambient  Temp.,  Tamb,  250.00  K;  Thermal  conductivity,  k,  2.50  watts/cm. K 
Disk  Gap  Factor,  g  =  2.0;  Number  of  numerical  intergration  steps,  n  «  300 

P  “  rO/rl,  lambda  -  i,  versus  radiation  fin  effectiveness 


\  p 

1.500 

3 . 000 

6 . 000 

10.000 

1  \ 

0.0001 

1.000 

1.000 

1 .000 

1 . 000 

0.0004 

0.999 

0.999 

0.999 

0.998 

0.0025 

0.994 

0.993 

0.991 

0.989 

0.0100 

0.97G 

0. 971' 

0.965 

0.958 

0.0225 

0.947 

0.938 

0.924 

0.911 

o 

o 

TT 

o 

o 

0.910 

0.895 

0.873 

0.853 

0.0900 

0.818 

0.793 

0.755 

0.724 

0. 1G00 

0.717 

0.685 

0.638 

0.599 

0.2500 

0.621 

0.585 

0.533 

0.493 

0.3G00 

0.533 

0.497 

0.446 

0.406 

0.4900 

0.458 

0.423 

0.374 

0.337 

0.6400 

0.394 

0.362 

0.317 

0.283 

0.8100 

0.341 

0.311  ' 

0.270 

0.239 

1.0000 

0.296 

0.270 

0.232 

0.205 
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p  -  rO/rl,  lambda  -  i,  versus  radiation  disk  OD  -cm. 


\  p 

1.500 

3 . 000 

6  .000 

10 . 000 

l  \ 

0.0001 

57.212 

143 . 036 

204  .352 

234.091 

0.0004 

14 .314 

35.791 

51.144 

58.600 

0.0025 

2.302 

5  .762 

8.246 

9.462 

0.0100 

0.586 

1.472 

2 .118 

2.442 

0.0225 

0.268 

0.678 

0.983 

1.141 

0.0400 

0.157 

0.399 

0.585 

0 .686 

0.0900 

0.078 

0.200 

0.300 

0.359 

0.1G00 

0.050 

0.130 

0.200 

0.244 

0.2500 

0.037 

0 . 098 

0.153 

0.190 

0.3600 

0.030 

0.080 

0.127 

0.160 

0.4900 

0.025 

0.069 

0 .111 

0.142 

0.6400 

0.023 

0 . 062 

0.101 

0.129 

0.8100 

0.021 

'0 . 057 

0.093 

0.121 

1.0000 

0.019 

0.053 

0.088 

0 .114 

P  «  rO/rl,  lambda  -  i,  versus  radiation  disk  ID  -  cm. 


\  P 
t  \ 

1.500 

3 .000 

6  .000 

10.000 

0.0001 

38.141 

47 .679 

34.059 

23.409 

0.0004 

9.542 

11.930 

8.524 

5.860 

0.0025 

1.535 

1.921 

1 . 374 

0.946 

0.0100 

0.391 

0.491 

0.353 

0.244 

0.0225 

0.179 

0.226 

0 . 164 

0 .114 

0.0400 

0.105 

0.133 

0.098 

0.069 

0.0900 

0.052 

.  0.067 

0.050 

0.036 

0.1600 

0.033 

0.043 

0.033 

0.024 

0.2500 

0.025 

0 . 033 

0 . 026 

0 . 019 

0.3600 

0.020 

0.027 

0.021 

0 .016 

0.4900 

0.017 

0.023 

0.019 

0  .  014 

0.6400 

0 . 015 

0.021 

0 . 017 

0 . 013 

0.8100 

0.014 

0.019 

0.016 

0 . 012 

1.0000 

0.013 

0 . 018 

0 . 015 

0  .  Oil 

P  “  rO/rl,  lambda  -  i,  versus  radiation  disk  thickness 


\  P 

1.500 

3  .000 

6 .000 

10  .  000 

1  \ 

0.0001 

9.627 

24  0 . 7  02 

767.655 

1174 .974 

0.0004 

0.151 

3.768 

12.021 

18.407 

0.0025 

0.001 

0 . 016 

0.050 

0 . 077 

0.0100 

0 .000 

0 .000 

0 .001 

0 .001 

0.0225 

0.000 

0 .000 

0.000 

0 .000 

0.0400 

0.000 

0.000 

0.000 

0 .000 

0.0900 

0.000 

0 .000 

0 .000 

0 . 000 

0.1600 

0.000 

0 . 000 

0 .000 

0 .000 

0.2500 

0.000 

0 . 000 

0 .000 

0.000 

0.3600 

0.000 

0 .000 

0  .  000 

0 .000 

0.4900 

0 . 000 

0 .000 

0.000 

0.000 

0.6400 

0.000 

0.000 

0 .000 

0 .000 

0.8100 

0.000 

,  0.000 

0.000 

0.000 

1.0000 

0.000 

0 .000 

0 . 000 

0.000 
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P  «  rO/ri,  lambda  -  i,  versus  Heat  Leak  Parameter,  Qs/QL 


\  p 

1.500 

3 .000 

6  .  000 

10 . 000 

X  \ 

0.0001 

1.5347 

'9 .5927 

19 .5797 

25  .6933 

0.0004 

0.0961 

-  0.6006 

1 .2264 

1.6100 

0.0025 

0 .0025 

0 . 0156 

0.0319 

0.0420 

0.0100 

0.0002 

0.0010 

0.0021 

0.0028 

0.0225 

0 . 0000 

0 . 0002 

0. 0005 

0.0006 

0. 0400 

0 . 0000 

0 .0001 

0 .0002 

0 .0002 

0.0900 

0.0000 

0.0000 

0.0000 

0 . 0001 

0.1600 

0.0000 

0.0000 

0 .0000 

0 . 0000 

0.2500 

0.0000 

0.0000 

0 .0000 

0 .0000 

0.3600 

0.0000 

0.0000 

0 .0000 

0 . 0000 

0.4900 

0.0000 

o.oooo 

0.0000 

0.0000 

0.6400 

0.-0000 

0.0000 

0.0000 

0 . 0000 

0.8100 

0.0000 

0.0000 

0.0000 

0 .0000 

1.0000 

0.0000 

0.0000 

0.0000 

0 . 0000 

p  -  rO/rl,  lambda  -  t,  versus  Rotor/Stator  Stack  Spacing 


\  p 
l  \ 

1.500 

3 .000 

6.000 

10.000 

0.0001 

28 .882 

722.105 

2302 . 966 

3524 .921 

0.0004 

0.452 

11.303 

36.064 

55.221 

0.0025 

0.002 

0.047 

0.150 

0.230 

0.0100 

0 . 000 

0 .001 

0.002 

0 . 004 

0.0225 

0.000 

0.000 

0 . 000 

0 . 000 

0.0400 

0.000 

0 .000 

0.000 

0 . 000 

0.0900 

0.000 

■  0.000 

0.000 

0 .000 

0.1600 

0 . 000 

0 .000 

0.000 

0 . 000 

0.2500 

0 . 000 

0 . 000 

0 . 000 

0  .  000 

0.3600 

0 . 000 

0 .000 

0 . 000 

0 . 000 

0.4900 

0.000 

0.000 

0 . 000 

0 . 000 

0.6400 

0.000 

0 .000 

0 . 000 

0 . 000 

0 .8100 

0.000 

0 . 000 

0 .000 

0  .  000 

1.0000 

0.000 

0.000 

0.000 

0 . 000 

P  »  rO/rl,  lambda  -  i,  versus  number  of  rotor  disks 


\  p 
l  \ 

0.0001 

1.500 

'  3 . 000 

6 .000 

10 . 000 

15 

2 

1 

1 

0,0004 

238 

24 

11 

8 

0.0025 

9231 

922 

4  13 

309 

0.0100 

144988 

14  4  35 

6  4  22 

4776 

0.0225 

712378 

70565 

31139 

22988 

0.0400 

2162582 

212836 

92975 

68  023 

0.0900 

9844698 

954507 

407409 

292143 

0.1600 

27295645 

2607049 

1087364 

764728 

0.2500 

57636064 

5432003 

2219775 

1535098 

0.3600 

102697103 

9574349 

3846821 

2623567 

0.4900 

163361787 

15098701 

5983707 

4  035828 

0.6400 

239837234 

22026497 

8634457 

5771795 

0.8100 

332183946 

30355956 

11798640 

7829863 

1.0000 

440301190 

40088133 

15473540 

10212739 
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RTC 
RTC Input 

Input  module  for  the  Rotary  Thermal  Coupler  design  program.  Note  that  the 
program  accepts  vector  strings  for  the  Radius  Ratio  and  the  Lamda  Input 
terms,  which  expedlts  parametric  studies.  The  terms  in  the  vector  strings 
are  to  be  seperated  by  commas.  Due  to  printout  limitations,  the  vector 
length  of  the  Radius  Ratio  parameter  is  limited  to  6  quantities. 

Enter  the  design  heat  load,  including  margin,  0.01  <  QL  <  20  watts;  1 

Enter  the  heat  source  temperature,  at  the  RTC  junction,  5K  <  TL  <  500K;  65 

Enter  the  heat  sink  temperature  at  the  RTC  interface,  4K  <  Tm  <  ( T 1— 1 ) ;  60 

Enter  the  ambient  environment  temperature,  TL  <  Tamb  <  500K  250 

Enter  the  disk  thermal  conductivity,  k  -  watts/cm. K;  2.5 

Enter  the  disk  gap  factor,  g,  where  gA=gap,  1  <  g  <  5;  2 

Enter  the  Module  Radius  Ratio  Vector,  1.5  <  P  =  DO/Di  <  10;  3.5 

Enter  the  Lamda  parameter  vector,  0.0001  <  i  <  1  -  0.0013,0.00135,0.0014 


RTCInitallzation 
RTCAna lysis 


i,  lamda; 

P*=Do/Dl ; 

Iterations 

0.001300 

3.500 

97 

0.001350 

3.500 

103 

0.001400 

3.500 

85 

Heat  Load  1.00  watts;  Source  temp.,  TL,  65.00  K;  Sink  temp.,  Tm,60.00K 

Ambient  Temp.,  Tamb,  250.00  K;  Thermal  conductivity,  k,  2.50  watts/cm. K 

Disk  Gap  Factor,  g  =  2.0;  Number  of  numerical  intergration  steps,  n  =  300 

P  =  rO/ri,  lambda  -  t,  versus  radiation  fin  effectiveness 

\  P  3.500 

l  \ 

0.00130  0.996 

0.00135  0.996 

0.00140  0.996 

P  =  rO/ri,  lambda  -  i,  versus  radiation  disk  OD  -  cm. 

\  P  3.500 

t  \ 

0.00130  12.271 

0.00135  11.818 

0.00140  11.398 
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\  p 
i  \ 
0.00130 
0.00135 
0.00140 


S  P 
t  \ 
0.00130 
0.00135 
0.00140 


\  P 
l  \ 
0.00130 
0.00135 
0.00140 


\  P 
l  \ 
0.00130 
0.00135 
0.00140 


\  P 
l  \ 
0.00130 
0.00135 
0.00140 


\  P 
l  \ 
0.00130 
0.00135 
0.00140 


p  *=  rO/rl,  lambda  -  i,  versus  radiation  disk  ID  -  cm. 

3.500 

3.506 

3.377 

3.257 

p  «=  rO/rl,1  lambda  -  i,  versus  radiation  disk  thickness  - 

3.500 

0.156 

0.140 

0.125  . 

P  -  rO/ri,  lambda  -  i,  versus  Heat  Leak  Parameter,  Qs/QL 
3.500 

0.0706 

0.0655 

0.0609 

P  •»  rO/ri,  lambda  -  i,  versus  Rotor/Stator  Stack  Spacing 

3.500 

0.469 

0.419 

0.376 

P  -  rO/ri,'  lambda  -  i,  versus  number  of  rotor  disks 

3.500 

197 

212 

228 

P  *=  rO/ri,  lambda  -  l,  versus  number  of  Stator  disks 

3.500 

•  198 
213 
229 


mm. 


mm 
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rtc 
RTC Input 

Input  module  for  the  Rotary  Thermal  Coupler  design  program.  Note  that  the 
program  accepts  vector  strings  for  the  Radius  Ratio  and  the  Lamda  Input 
terms,  which  expedlts  parametric  studies.  The  terms  in  the  vector  strings 
are  to  be  seperated  by  spacebar  blanks.  Due  to  printout  limitations,  the  vector 
length  of  the  Radius  Ratio  parameter  is  limited  to  6  quantities. 

Enter  the  design  heat  load,  including  margin,  0.01  <  QL  <  20  watts;  ,35 

Enter  the  heat  source  temperature,  at  the  RTC  Junction,  5K  <  TL  <  500K;  40 

Enter  the  heat  sink  temperature  at  the  RTC  interface,  4K  <  Tm  <  ( Tl— 1 ) ;  35 

Enter  the  ambient  environment  temperature,  TL  <  Tamb  <  500K  120 

Enter  the  disk  thermal  conductivity,  k  -  watts/cm. K;  2.5 

Enter  the  disk  gap  factor,  g,  where  gA>=gap,  1  <  g  <  5 ;  2 

Enter  the  Module  Radius  Ratio  Vector,  1.5  <  P  “  D0/D1  <  10;  3, 3. 5, 4, 4. 5 

Enter  the  Lamda  parameter  vector,  0.0001  <  i  <  1  -  . 00035 ,. 00055 ,. 00075 ,. 001 ,. 0015 ,. 002 

RTCInltallzat ion 
RTCAna lysis 

Heat  Load  0.35  watts;  Source  temp.,  TL,  40.00  K;  Sink  temp.,  Tm,35.00K 
Ambient  Temp.,  Tamb,  120.00  K;  Thermal  conductivity,  k,  2.50  watts/cm. K 
Disk  Gap  Factor,  g  -  2.0;  Number  of  numerical  intergratlon  steps,  n  “  300 

P  >•  rO/ri,  lambda  -  t,  versus  radiation  fin  effectiveness 


\  p 

3 .000 

3.500 

4.000 

4.500 

l  \ 

0 .00035 

0.999 

0.999 

0.999 

0.999 

0 . 00055 

0.998 

0.998 

0.998 

0.998 

0.00075 

0.998 

0.998 

0.998 

0.998 

0.00100 

0.997 

0.997 

0.997 

0.997 

0.00150 

0.996 

0.996 

0.995 

0.995 

0.00200 

0.994 

0.994 

0.994 

0.994 

p  -  rO/ri 

,  lambda 

-  i,  versus 

radiation  disk 

\  P 

3.000 

3.500 

O 

O 

o 

4.500 

1  \ 

0.00035 

15.399 

17.111 

18.481 

19.602 

0 . 00055 

9.805 

10.896 

11.768 

12.482 

0.00075 

7.195 

7.995 

8.635 

9.160 

0.00100 

5.400 

6 . 001 

6.482 

6  .875 

0.00150 

3.605 

4  .  006 

4  .328 

4.591 

0.00200 

2.708 

3 . 009 

3.251 

3.449 
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P  -  rO/ri,  lambda  -  i ,  versus  radiation  disk  ID  -  cm. 


\  p 

3.000 

3.500 

4  .000 

4.500 

l  \ 

0.00035 

5.133 

4.889 

4.620 

4.356 

0.00055 

3.268 

?  .  113 

2.942 

2 .774 

0.00075 

2.398 

2.284 

2.159 

2  .  035 

0.00100 

1.800 

1.714 

1.620 

1.528 

0.00150 

1.202 

1.145 

1.082 

1.020 

0.00200 

0.903 

0.860 

0.813 

0.766 

p  rO/ri, 

lambda 

-  i,  versus 

radiation  disk  thickness  - 

\  P 

3 . 000 

3.500 

4  .000 

4.500 

1  \ 

0.00035 

0.186 

0.263 

0.339 

0.410 

0.00055 

0.048 

0.068 

0.08  7 

0.106 

0.00075 

0.019 

0.027 

0 . 035 

0.  042 

0.00100 

0.008 

0.011 

0 . 015 

0  .  018 

0.00150 

0.002 

0.003 

0 . 004 

0  .  005 

0.00200 

0.001 

0. 001 

0 . 002 

0.002 

p  =  rO/ri, 

lambda 

-  i  ,  versus 

Heat  Leak  Parameter,  Qs/QL 

\  P 

3 . 000 

3.500 

4  .000 

4.500 

l  \ 

0.00035 

0.0168  0 

.  0207 

0.0242 

0 . 0272 

0.00055 

0.0068  0 

.  0084 

0.0098 

0.0110 

0.00075 

0.0037  0 

.  0045 

0.0053 

0.0059 

0.00100 

0.0021  0 

.  0026 

0.0030 

0 . 0033 

0.00150 

0.0009  0 

.,0011 

0 .0013 

0 . 0015 

0.00200 

0.0005  0 

.  0006 

0 . 0007 

0 . 0008 

p  =  rO/rl, 

lambda 

-  i ,  versus 

Rotor/Stator  Stack  Spacing 

\  P 

3 . 000 

3.500 

4.000 

4.500 

i  \ 

0.00035 

0.557 

0.790 

1 . 016 

1.229 

0.00055 

0.144 

0.204 

0.262 

0.317 

0.00075 

0.05  7 

0.080 

O 

F- 1 
O 

0.125 

0.00100 

0.024 

0 . 034 

0.044 

0.053 

0.00150 

0.007 

0.010 

0.013 

0  .  016 

0.00200 

0.003 

0 . 004 

0.006 

0 .007 

P  =  rO/ri, 

lambda 

-  i,  versus 

number  of  rotor  disks 

\  P 

3 .000 

3.500 

4  .  000 

4.500 

l  \ 

0 .00035 

208 

163 

137 

120 

0.00055 

512 

4  01 

337 

296 

0.00075 

951 

746 

626 

54  9 

0.00100 

1689 

1324 

1112 

975 

0.00150 

3794 

2974 

24  9  8 

2189 

0.00200 

6736 

5279 

4  4  33 

3885 

mm. 


mm 


Electro  Thermo  Associates 
60 


•  Phase  I  Final  Report 

Rotary  Cryogenic  Thermal  Coupling 


RTc 

RTC Input 

Input  module  for  the  Rotary  Thermal  Coupler  design  program.  Note  that  the 
program  accepts  vector  strings  for  the  Radius  Ratio  and  the  Lamda  input 
terms,  which  expedits  parametric  studies.  The  terms  in  the  vector  strings 
are  to  be  seperated  by  spacebar  blanks.  Due  to  printout  limitations,  the  vector 
length  of  the  Radius  Ratio  parameter  is  limited  to  6  quantities. 

Enter  the  design  heat  load,  including  margin,  0.01  <  QL  <  20  watts;  .15 

Enter  the  heat  source  temperature,  at  the  RTC  junction,  5K  <  TL  <  500K;  12 

Enter  the  heat  sink  temperature  at  the  RTC  interface,  4K  <  Tm  <  (Tl-1);  10 

Enter  the  ambient  environment  temperature,  TL  <  Tamb  <  500K  60 

Enter  the  disk  thermal  conductivity,  k  -  watts/cm. K;  1 

Enter  the  disk  gap  factor,  g,  where  gA=gap,  1  <  g  <  5;  2 

Enter  the  Module  Radius  Ratio  Vector,  1.5  <  p  =  D0/D1  <10;  345 

Enter  the  Lamda  parameter  vector,  0.0001  <  *  <  1  -  . 0002 ,. 0003 ,. 0004 ,. 0005 

RTCInltalizatlon 

RTCAnalysis 

Heat  Load  0.15  watts;  Source  temp.,  TL,  12.00  K;  Sink  temp.,  Tm,10.00K 
Ambient  Temp.,  Tamb,  60.00  K;  Thermal  conductivity,  k,  1.00  watts/cm. K 
Disk  Gap  Factor,  g  =  2.0;  Number  of  numerical  intergration  steps,  n  =  300 

P  =  rO/ri,  lambda  -  i,  versus  radiation  fin  effectiveness 


\  p 

3 . 000 

4  .000 

5  .000 

l  \ 

0.00020 

0.999 

0.999 

0.999 

0.00030 

0.999 

0.999 

0.999 

0.00040 

0.999 

0 .999 

0.999 

0.00050 

0.999 

0.998  1 

0.998 

P  =  rO/ri,  lambda  -  i,  versus  radiation  disk  OD  -  cm. 


\  P 

3 . 000 

4  .000 

5  .000 

l  \ 

0.00020 

76.893 

92.277  . 

102.535 

0.00030 

51.276 

61.537 

68.379 

0.00040 

38.468 

46 . 167 

51.301 

0.00050 

30.783 

36.945 

41.055 
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p  =  rO/ri 

,  lambda 

-  x ,  versus 

radiation  disk  ID  -  cm. 

\  p 

3 . 000 

4  .000 

5 .000 

t  \ 

0.00020 

25.631 

23.069 

20.507 

0.00030 

17 . 092 

15 .384 

13 .676 

0.00040 

12.823 

11.542 

10.260 

0.00050 

10.261 

9.236 

8.211 

P  =  rO/ri,  lambda 

l 

-  i,  versus 

radiation  disk  thickness  - 

\  P 

3 .000 

4  .000 

5 .000 

X  \ 

0.00020 

0.547 

0.997 

1.401 

0.00030 

0.162 

0.296 

0.415 

0.00040 

0.068 

0.125 

0.175 

0.00050 

0.035 

0.064 

0.090 

P  -  rO/rl 

,  lambda 

-  i,  versus 

Heat  Leak  Parameter,  Qs/QL 

\  P 

3 .000 

4  .000 

5 .000 

l  \ 

0.00020 

0.0615 

0.0885 

0.1093 

0.00030 

0.0273 

0.0394 

0.0486 

0.00040 

0.0154 

0.0222 

0 . 0274 

0.00050 

0 . 0099 

0.0142 

0.0175 

P  *  rO/ri 

,  lambda 

-  x,  versus 

Rotor/Stator  Stack  Spacing 

\  P 

3 . 000 

4  .  000 

5  .  000 

i  \ 

0.00020 

1.641 

2.992 

4.203 

0.00030 

0.487 

0 .887 

1.246 

0.00040 

0.205 

0.374 

0 .526 

0.00050 

0.105 

p  .192 

0.270 

P  «  rO/rl 

,  lambda 

-  x,  versus 

number  of  rotor  disks 

\  P 

3 . 000 

4.000 

5  .  000 

X  \ 

0.00020 

352 

232 

183 

0.00030 

791 

521 

412 

0. 00040 

1405 

925 

732 

0.00050 

2195 

1445 

1143 

P  =  rO/ri 

,  lambda 

-  x,  versus 

number  of  Stator  disks 

\  P 

3.000 

4  .  000 

5 .000 

l  \ 

0.00020 

353 

233 

18  4 

0.00030 

792 

,  522 

413 

0.00040 

1406 

926 

733 

0.00050 

2196 

1446 

114  4 

nun . 


nun 
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